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Abstract 

We examine the use of adiabatic quantum algorithms to solve struc- 
tured, or nested, search problems. We construct suitable time depen- 
dent Hamiltonians and derive the computation times for a general class 
of nested searches involving n qubits. As expected, we find that as ad- 
ditional structure is included, the Hamiltonians become more local and 
the computation times decrease. 



1 Introduction 



Adiabatic Quantum Computation (AQC) is a relatively new paradigm in the 
field of quantum computing. Whereas in the standard model of QC (SQC), 
an algorithm is defined as a sequence of discrete unitary transformations, 
AQC (see, e.g., ||, [2], ||, ||, ||, |6|) considers the continuous time evolution of 
the quantum system, described by the Schrodinger equation 
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with a time dependent Hamiltonian H{t). A computational problem which 
has been studied extensively in both SQC and AQC is the unstructured 
search. In its simplest form the unstructured search corresponds to find- 
ing a single marked state (needle) in a completely unstructured database 
(haystack) of iV states. In the case that the physical system consists of n 
two state particles (qubits), N = 2 n corresponds to the size of the complete 
Hilbert space for the system. Classically a random search of N objects re- 
quires on average N operations to pick out the marked object. It was shown 
by BBBV @ that for SQC, the lower bound on the corresponding number 
of steps is of order y/~N. Grover was the first to construct a specific QC 
algorithm that achieved this lower bound M. 

In the context of AQC, the physical quantity that seems to most closely 
correspond to the number of steps of a SQC algorithm is the time required to 
do the computation. It is therefore not surprising that under fairly standard 
assumptions, it can be shown that the time required for the unstructured 
search in AQC increases as []. 

In many physical systems, it is possible to use additional structure in 
order to shorten the search. For example, consider a database of items with 
n distinguishable bits. A maximally structured classical search would check 
one bit at a time in order to find the correct marked state. Such a search 
would require O(n) operations. The key is that, by using the existing struc- 
ture, one is able in principle to search n two-dimensional objects, instead of 
a random search among N = 2 n objects. 

The above example is a special case of a general class of nested searches, 
which have recently been analyzed || in the context of SQC. The purpose of 
the present paper is to study structured adiabatic quantum search. We find 
that additional structure shortens the running time and makes the Hamil- 
tonian more local (concerning the interactions of the system, as is explained 
below). Our results for the running time are consistent with the conjec- 
ture (cf. ||) that (adiabatic) quantum computation improves the running 
time of the corresponding classical algorithm by a square root; due to the 
fact that quantum computation manipulates probability amplitudes, whose 
squares give the corresponding probabilities. 

We start with a short description of AQC in general. In Section 3 we 
show that the unstructured search studied by Grover leads to a spatially 
non-local Hamiltonian when considered in the conventional model for the 

1 In a dynamical quantum system one can always shorten the time by increasing the 
energy. The consequences for AQC of increasing the energy temporarily have recently 
been analyzed in 
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implementation of quantum algorithms, namely two-level quantum systems 
(qubits), such as spin- 1/2 particles. That is, the unstructured adiabatic 
search requires an ra-body interaction for an n-qubit system. In Section 4, we 
define the notion of a structured adiabatic quantum search Hamiltonian and 
show that it is in general local, in the sense referred to above: it couples only 
a subset of the total number of qubits. We then start our detailed analysis 
of nested searches with the most intuitive case, the maximally structured 
search, and consider its adiabatic evolution and running time, which is shown 
to scale as T = 0{y/n). Finally, we study the general case and obtain its 
running time, comparing different ways of structuring the system. 



2 Adiabatic Quantum Computation 

The method of AQC is based on the use of the adiabatic theorem for con- 
sidering quantum computing as continuous time evolution from some easily 
prepared initial state \ipi) to a final state \ipf) that encodes the solution to 
the computational problem. Specifically, the adiabatic theorem (see, e.g. 



[ 10 1 ) states the following: a system which is described by a time-dependent 
Hamiltonian H (t) will stay close to the instantaneous ground state of H(t) 
provided that the time evolution is slow enough. In quantitative terms, after 
time T 

\(E ;T\y(T))\ 2 >l-e 2 , (2) 
where Eq is the lowest energy eigenvalue and e <C 1, provided that: 

K^fl^U. (3) 



mm 



where 



uw = min \E x {t) - E (t)} (4) 
0<t<T 



is the minimum gap between the lowest two energy eigenvalues Eq and E\. 
Thus, if (p|) is satisfied, a measurement after time T gives as result the 
solution with almost certainty, i.e. with probability ~ 1 — e 2 . 

The time-dependent Hamiltonian that is traditionally used in AQC is 
constructed as the linear combination 

H(t)=f(t)H i + g(t)H f (5) 

where \ipi), \ipf) are the ground states of Hi, Hf, respectively and f(t) and 
g(t) are usually considered to be monotonic functions such that /(0) = 1, 
f(T) = 0, g(0) = and g(T) = 1, where T is the total computation time. 
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It will be useful for what follows to note that the adiabatic theorem (||) 
can be modified [^] for the case of an m-fold degenerate first excited state. 
In this case, it reads 

m I / 771 I dH I 771 \ 1 2 

V '^'T 1 0)t ' <e 2 - (6) 

i=l mm 

3 Unstructured search in AQC 

The goal of the search algorithm is to find a marked object \m) in an un- 
structured database of size N in as few steps as possible. Classically, one 
has to look at O(N) objects to find the marked one. One advantage of the 
quantum search is that it manipulates the amplitude of the marked state, 
leading to a quadratic amplification of the probability. The initial state is 
the superposition of all states with equal weight, 

N 1 

* ) = i<i/(o)> = ]r— 1<>. (7) 

i=i viv 



A natural choice for a Hamiltonian that has this state as its ground state is: 

fli = l-|*o>(*o|. (8) 

Correspondingly one must choose (cf. Section 4.1), 

H f = l-\m)(m\, (9) 

which has the ground state \m). Hf is in a sense the adiabatic analogue of 
the "oracle" I — 2\m){m\ used in the Grover search algorithm. It projects 
out the marked state \m). Thus, in the framework of AQC, the initial state 
|\&o) evolves into \m) in time T, given that the Hamiltonian (|J) varies slowly 
enough with time. In the simplest case f(t) = 1 — t/T and g(t) = t/T, the 
minimum running time increases with N at the same rate as for the classical 
search: T = 0(N). However, one can improve this by choosing f(t) and g(t) 
to vary most rapidly when the gap uj is the largest. For example, one can 
use the approach of in which f(t) = 1 — s(t), g(t) = s(t) with s(0) = 
and s(T) = 1. This choice yields an adiabaticity condition of the form: 



K£ i;g |^ft|go; g )l 
w 2 (s) 



ds 



<e. (10) 



One is then free to choose s(t) so that the bound in (10) is saturated for all 



t. This yields a running time T = O(yN). Recently, it was shown [11] that 
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with a more general choice of the functions / and g in the construction of 
H(t) the running time can - in principle - be made independent of N, but 
at the cost of requiring a large amount of energy to be temporarily injected 
into the system. 

We are now lead to the question of how to implement a generic AQC 
algorithm in a real quantum computer, and what kind of physical systems 
might be suitable. The conventional scheme for the implementation of quan- 
tum computation considers a system of uncorrelated 2-level quantum sys- 
tems (qubits), e.g. spin 1/2-particles. The 1-qubit space TLi is spanned by 
|0) = (q) and |1) = The n-qubit space is given by the tensor product 

n = Hi ®n 2 ® ... ® n n (11) 

with the basis {\ati, a n ) = \ati)® ...(8) |a n )|aj € {0, 1}}. The marked state 
is \m) = \z\...z n ) = \zi) ••• <8> \z n ) with Z( £ {0, 1}. In this framework, the 
final Hamiltonian 

H f = l- \m)(m\ (12) 
can be written in terms of Pauli spin matrices 

H f = i-n^i-HM')) as) 

n i 

= E^n^- ( i4 ) 

In this form it is clear that the Hamiltonian contains a product of n spin 
matrices, which implies interactions among all qubits. This non-locality may 
put constraints on possible implementations of such algorithms, at least for 
spin systems^. 

In the next section, we describe how, by using additional structure, the 
search can be modified to reduce the spatial non-locality in the Hamiltonian, 
while keeping the initial and final ground state the same. This structured 
search also provides a shorter computation time than the unstructured one. 



4 Structured adiabatic search 

The adiabatic theorem requires the evolution of the quantum system to be 
slow enough to keep the system in its ground state for all times. It seems 

2 If one considers AQC only as a possible simulation of SQC, the non-locality is irrele- 
vant. 
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plausible that it is "easier" for the system to stay in its ground state if there 
are no long range interactions, but only 2 or 3-qubit interactions. In order 
to find a Hamiltonian that realizes this idea, it is useful to consider the 
problem in the language of decision clauses etc.(cf. ||): 

4.1 Satisfiability problems, clauses and Hamiltonians 

Search algorithms belong to the class of satisfiability problems. A satisfia- 
bility problem is a combination of decision clauses Cj which have an element 
of {0, 1} (true or false) as output which depends on the values of some set 
of the qubits of the system. A formula is an rz-bit instance of satisfiability, 



i.e. in general m clauses acting on n bits. The action of a clause or formula 
on an assignment can be translated into the action of a Hamiltonian on the 
states of a physical system (up to an overall factor): Formula ( |15| ) is mapped 
into 



where each Hi is constructed from the clause and acts only on the bits of 
this clause. For example, if C{ is a two-bit clause, Hi contains only operators 
that act on these two bits, it can have at most a two-body interaction term. 
The action of Hi on the bits absent in Cj is given by the identity operator. 
The initial state is the ground state of H(0), and the ground state of H(t = 
T) satisfies all the clauses in the formula. 

In this language, the unstructured search algorithm is a single clause that 
acts on n qubits and has a unique (but unknown) satisfying assignment, \m). 

It should be noted that the form of Hf given in (||) is the most general 
one available for an oracle search. By definition, an oracle clause acts as 



F = d A C 2 A ... A C, 



(15) 



H(t) = H 1 (t) + H 2 (t) + ... + H n (t) 



(16) 








(17) 
(18) 




where f(T) should be the same for all i / m. Therefore, Hf is diagonal, 
and it reads 



Hf = /(T)£|i)<i 
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4.2 Maximally structured search 



The unstructured search, which is non-local in space, corresponds to a single 
n-bit clause. In order to find a local Hamiltonian, it is therefore the most 
intuitive alternative to consider a formula of n 1-bit clauses, 



F = d AC 2 A... AC n , 



(22) 



where C% is satisfied if and only if the i-th bit has the required value Z{ 
for the marked state \m) = \z\Z2---z n ). We call this "maximally structured" 
because it corresponds to the maximal splitting of the system: The adiabatic 
Hamiltonian which has as ground state the unique satisfying assignment of 
F is a sum of 1-bit versions of the oracle Hamiltonian, 



H 



sir 



i=i 



(23) 



with 



hi 



(l-\Zi)(Zi\)®l 



1® 



(1 _ (-1)^(0) , 



(24) 



which act as oracle on the i-th bit and as identity on the others. The 
application of Hj tr on a general state \a\a2---a n ) gives the number of qubits 
with oti ^ Zi. This illustrates the crucial difference between this procedure 
and the unstructured search. The latter leads to the same answer, 1, for all 
unsatisfying assignments, whereas in the present case the query gives the 
number of unsatisfied clauses. Thus, one gets additional information. 

As before, the initial state |^o) is the ground state of Hi = 1 — l^oX^ol) 
which can be written as 



H 



sir 



After some algebra (cf. |l 

H str (t) 
are found to be 



the eigenvalues of the Hamiltonian 

I ii n 

i=i i=i 



E 



str 



n 



(f + g)-mJp + g 



(25) 

(26) 
(27) 
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where m = J2 m ii with m i = ±1/2 representing the spin up and down of 
the i-th particle, corresponding to ai = 0, 1. The ground state is obtained 
for rrii = 1/2 Vi, i.e. m = n/2. The first excited state (m = n/2 — 1) is 
n-fold degenerate. The crucial condition for the adiabatic theorem to hold 
is a non-vanishing gap between the ground state and the first excited state, 
which is fulfilled: 



W Ei?f''-Bf = ^f+ffV0. (28) 

Note that this result corresponds to Eq.(14) of [ O ] for 2™ = N = 2, i.e. for 
n = 1, the single qubit adiabatic search. The running time is determined by 
the adiabatic theorem, which in the case of a degenerate first excited state 
leads to the condition (|6|). The transition probability to each of the first 
excited states is proportional to the square of 

{E{tr] dH = _(-!)« fg-gf 

* N 1 1 dt 1 1 2 y/f 2 +g 2 



so the total transition probability is n times Eq.(29) squared (cf. Eq.(15) 



in [11] for N = 2). For an estimate of the running time, we consider the 
simplest case in which / and g are linear in s, f(s) = l — s(t) and g(s) = s(t), 
which, after optimizing s(t) as in the previous section, results in (cf.]|l2]l) 



T = vVe = VlogN/e , (30) 

polynomial in n. It should be noted that a similar polynomial time algorithm 
has been studied in ||]. 

4.3 General case of structured search 

The transformation of the n-bit clause into n 1-bit clauses is not the only 
possibility of introducing structure to the search. For example, one could 
also consider clauses acting on 2 qubits, 

n— l)n 

(31) 

corresponding to the Hamiltonian 

n/2-1 

H str = ^ /l[2i+l][2(i+l)] » ( 32 ) 
i=0 

where 

h[2i+l][2(i+l)] = 1<8>...<8>1® (1- \Z2i+lZ2(i +1 ))(Z2i + iZ2( i+ i)\)®l®...®l (33) 
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m 


n/m 


eT 


a 


P 


1 


6 


7.94 


0.9962 


oo 


2 


3 


3.74 


0.9518 


3.8074 


3 


2 


3.00 


0.8842 


2.0000 


6 


1 


2.45 


0.7211 


1.0000 



Table 1: Comparison of times for various splittings for n = 6. 



m 


n/m 


eT 


a 


P 


1 


30 


32768.00 


1.0000 


oo 


2 


15 


256.00 


1.0000 


16.0000 


3 


10 


55.40 


0.9999 


7.3084 


5 


6 


17.75 


0.9973 


3.5743 


6 


5 


13.64 


0.9940 


2.9165 


10 


3 


8.37 


0.9695 


1.8451 


15 


2 


6.71 


0.9297 


1.4057 


30 


1 


5.48 


0.8307 


1.0000 



Table 2: Comparison of times for various splittings for n = 30. 

acts on the two neighboring qubits as oracle and as identity on the others. 
A more interesting case with respect to realization would be 

n-l 

Hf = £ h i[i+1] , (34) 
i=i 

with overlap of the interacting qubit-pairs (see the discussion in [|l^] ) . 

In general, any splitting of the N = 2 n -dimensional Hilbert space into m 
smaller spaces of dimensions iVj = 2 ni with n % = n is possible. In order to 
compare the effects of different ways of splitting, we calculate the running 



time for several cases numerically (cf. [11 1 for the method). 

In order to investigate the scaling of T with the splitting, we define the 
two coefficients a, P implicitly by: 



eT = (v^vV/™)" 

eT = (vW 3 • (35) 

We find (see Tables 1 and 2) that for no splitting (m = 1), a — ► 1 for 
large n, which means the running time is T = 0(y/~N), as is known from the 
unstructured search. 
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For m = 2, i.e. n = n\ + U2, the time is shortened, and the optimal value 
in this case is achieved for m = 112- 

For m > 2, the running time is even shorter, and it scales with the root 
of the dimension of the largest Hilbert space of the splitting, so it is again 
optimal for equal values of n$ for all i. This can also be seen from the 
expression for the running time, 



1 r l 



with 



e Jo 



ds\fg-gf\ t 



^ Nj - 1 1 



\f-9? + W J9 



(36) 



(37) 



the gap for the i-th subsystem. For splitting into m parts of equal size, 
N . = 2 n/m Vij and Eq.@) simplifies to 



1 vw- 

T = -yjm — 

e Ni 



ds 



\fg-gf\ 



[(/ - 9? + 4-fg?/ 2 



(38) 



(which reduces to Eq.(20) in [11] for m = 1) where the integral scales with 
Ni. For n/m 3> 1, we find a ~ 1, and 



eT 



For maximal splitting, m = n and n 
results in T = 0(y/n), as shown above. 



m \Jl n l m . (39) 
1 Vi, leading to (3 = 1, which 



5 Conclusions 

We have shown how additional structure affects both the running time and 
locality of the Hamiltonian required for an adiabatic quantum search algo- 
rithm. As expected, the more structure, the more local the Hamiltonian 
and the shorter the running time. In fact, the expressions we have obtained 
for the running time of the general structured search suggest strongly that 
adiabatic quantum computation consistently improves the running time for 
the corresponding classical algorithm by a square root. This supports the 
conjecture (cf. ||) that the speed-up achieved by quantum computation can 
be directly attributed to the fact that the quantum computation algorithm 
manipulates probability amplitudes, whose squares give the corresponding 
probabilities. It would be interesting to see whether this conjecture can be 



10 



proven in more general contexts. Another interesting topic for future work 
concerns what happens in the more general case that the pairwise interac- 
tions overlap. These and other related issues will be addressed elsewhere. 
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